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ABSTRACT 


A  multiple-scattering  laser  propagation  model  is  developed  in 
the  form  of  three  simultaneous  partial  differential  equations  for  the 
reduced  field  amplitude,  the  forward-  and  the  backscattered  irradi- 
ances.  The  scattering  and  diffusion  coefficients  are  defined  in  terms 
of  Mie  functions,  and  reduction  factors  are  derived  to  take  into 
account  detection  at  limited  fields  of  view.  Solutions  for  transmitted 
beam  profiles  and  integrated  backscatter  are  compared  with  measurements 
performed  on  laboratory*generated  water  droplet  clouds  at  0.63,  1.06 
and  10.6  pm.  The  agreement  is  very  good  for  all  conditions  investi¬ 
gated,  ranging  from  negligible  to  important  multiple-scattering 
contributions.  J  1 


RESUME 


Un  modele  de  propagation  laser  en  presence  de  diffusions  multi¬ 
ples  par  des  particules  en  suspension  est  mis  au  point  sous  la  forme  de 
trois  Equations  aux  derivees  partielles  pour  l'amplitude  reduite  du 
champ,  1' intensity  lumineuse  dif fusee  vers  l'avant  et  1* intensity 
lumineuse  diffusee  vers  l’arriSre.  Les  coefficients  d'attenuation  et 
de  diffusion  sont  dSfinis  au  moyen  de  fonctions  de  Mie,  et  des  facteurs 
de  reduction  sont  obtenus  pour  tenir  compte  du  champ  de  vue  limits  des 
dStecteurs.  Les  solutions  du  profil  des  faisceaux  transmis  et  de  la 
rfctrodiffusion  sont  comparSes  t  des  mesures  faites  en  laboratoire  3 
0.63,  1.06  et  10.6  pm  sur  des  nuages  de  gouttelettes  d'eau.  L’accord 
est  tres  bon  pour  tous  les  cas  etudies,  les  contributions  des  diffu¬ 
sions  multiples  variant  de  nSgligeables  3  importantes. 
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1.0  INTRODUCTION 

In  recent  years,  many  electro-optical  systems  have  been  designed 
to  enhance  the  effectiveness  of  various  military  weapons.  In  clear 
weather,  these  systems  are  very  predictive  and  accurate.  However,  poor 
weather  conditions  and/or  artificial  obscurants  can  severely  degrade 
their  performance.  Thus,  a  requirement  exists  for  understanding  and 
predicting  the  propagation  of  electromagnetic  waves  under  low- 
visibility  conditions  created  mostly  by  the  presence  of  natural  or  man¬ 
made  aerosols.  Also,  theoretical  developments  are  needed  for  the 
inverse  problem  of  determining  the  atmospheric  optical  and  Infrared 
parameters  from  single-ended  lidar  backscatter  measurements. 

Propagation  in  thin  or  tenuous  clouds  is  well  understood  and 
documented  (e.g.  Refs.  1~2).  It  is  governed  by  simple  and  exact 
single-scattering  equations  which  explain  numerous  atmospheric  phenom¬ 
ena.  However,  in  dense  clouds,  the  single-scattering  approximation 
fails.  The  optical  rays  undergo  many  scattering  events  before  escaping 
the  medium,  reaching  a  target,  or  being  detected.  These  multiple  scat¬ 
terings  have  nonnegligible  and  often  dominant  effects  on  the  transmit¬ 
ted  and  backscattered  waves.  Analytic  methods  (e.g.  Refs.  3-5), 
transport  methods  (e.g.  Refs.  6-12),  and  Monte  Carlo  simulations  (e.g. 
Refs.  13-14)  have  been  used  to  study  this  problem.  They  all  constitute 
very  powerful  tools,  but  the  mathematics  remains  complex  and  no  general 
solutions  have  been  obtained  yet. 

In  a  preceding  report  (Ref.  15),  we  set  the  foundations  of  a 


multiple-scattering  irradiance  propagation  model  based  on  a  stochastic 
representation  of  the  refractive-index  field  created  by  the  presence  of 
randomly  distributed  suspended  particles.  In  the  present  report,  we 
use  the  propagation  equations  derived  in  Ref.  15  and  propose  a 
heuristic  but  convenient  determination  of  the  propagation  coefficients. 
Also,  we  introduce  f ield-of-view  functions  to  take  into  account 
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practical  detector  geometries.  The  results  are  in  good  agreement  with 
forward-  and  backscattering  measurements  in  laboratory-generated  water 
droplet  clouds. 

The  propagation  equations  are  given  in  Chapter  2.0  and  the 
coefficients,  in  Chapter  3.0.  The  solutions  are  derived  in  Chapter  4.0 
and  the  field-of-view  functions,  in  Chapter  5.0.  Chapter  6.0  compares 
the  model  predictions  with  laboratory  simulation  data. 

This  work  was  performed  at  DREV  between  April  1985  and  January 
1986  under  PCN  33B13,  Multiple-Scattering  Effects  in  Atmospheric 
Aerosols. 


2.0  PROPAGATION  MODEL 


The  model  stems  from  the  wave  equation  of  electromagnetic  propa¬ 
gation.  In  the  version  discussed  in  this  report,  we  do  not  consider 
the  polarization  effects  and  we  restrict  the  analysis  to  one  temporal 
Fourier  component  of  the  wave.  We  thus  proceed  from  the  temporal 
Fourier  transform  of  the  scalar  wave  equation,  i.e. 

V2E  -  E  -  S,  [!] 

C2 

where  E  is  the  Fourier  transform  of  the  scalar  electric  field,  V2  is 
the  three-dimensional  Laplacian  operator,  n  is  the  refractive  index,  w 
is  the  angular  frequency  of  the  wave,  c  is  the  speed  of  light  in  free 
space,  and  the  right-hand  side  function  S  is  a  source-sink  term. 

Equation  1  is  treated  as  a  stochastic  equation  since  the  random 
spatial,  temporal  and  size  distributions  of  the  suspended  particles 
constitute  a  space-time  random  field  for  the  complex  refractive  index 
n.  The  statistical  properties  of  the  function  n  can  be  derived  from 
the  number  density,  shape,  size  distribution  and  complex  refractive 
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index  of  the  particles,  and  the  refractive  index  of  the  surrounding 
medium.  The  model  actually  postulates  two  equations  similar  to  eq.  1 
to  separate  the  electric  field  into  a  forward  and  a  backward  wave.  The 
two  waves  Interact  through  their  respective  source-sink  term  S  which 
accounts  for  the  backscatterings  from  the  particles.  The  formulation 
of  the  source-sink  terms  constitutes  a  basic  feature  of  the  model  and 
is  discussed  at  length  in  Ref.  15.  Here,  we  only  reproduce,  with  a 
slight  change  in  notation,  the  resulting  equations  for  the  forward  and 
backward  field  components  E*  and  E“,  i.e. 

V2E+  -  E+  -  ik  e21kZ  E-[v2*“  -  V2*+],  [2  ] 

c2 

V2E“  -  E"  -  ik  e  -2ikz  E+fv2^  -  V2 4>“  ] ,  [3] 

c2 

where  k  ■  n  w/c  is  the  wave  number  in  the  surrounding  medium  of  index 
o 

nQ,  z  is  the  coordinate  along  the  main  direction  of  propagation,  is 
the  random  geometrical  phase  (expressed  in  units  of  length)  of  the  wave 
backscattered  from  the  forward-propagating  component  E+,  and  4>~  is  the 
random  phase  of  the  wave  backscattered  from  the  backward-propagating 
component  E~. 

The  mechanism  of  separation  of  the  electromagnetic  wave  into  a 
forward  and  a  backward  component  modeled  by  eqs.  2  and  3  is  arbitrary. 
The  original  Maxwell's  equations  contain  no  such  distinction.  However, 
E+  and  E~  constitute  a  legitimate  representation  of  the  total  field. 

One  can  always  rewrite  a  solution  in  multiple  parts  by  imposing  addi¬ 
tional  conditions.  The  condition  here  Is  the  source-sink  function  or 
the  right-hand  side  term  in  eqs.  2  and  3.  This  separation  Is  very 
helpful  in  practice  as  it  leads  to  solutions  that  are  experimentally 
distinguishable. 
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The  stochastic  equations  1  and  2  are  used  to  obtain,  by  standard 
and  straightforward  operations,  the  deterministic  equations  for  the 
first-  and  second-order  statistical  moments  of  E+  and  E~,  l.e.  the 
average  field  amplitude  and  the  average  irradiance  of  the  forward  and 
backward  waves.  This  produces  an  unclosed  hierarchy  of  equations  since 
there  are  more  unknowns  than  equations.  The  essence  of  the  model 
consists  In  the  derivation  of  closure  relations  to  link  the  higher 
order  unknown  moments  to  the  lower  order  average  field  amplitude  and 
irradiance  functions.  This  Is  accomplished  through  the  hypotheses  of 
quasl-homogenelty  and  quasi-isotropy  of  the  random  amplitude  and  phase 
functions  and  of  weak  but  nonnegliglble  correlation  between  the  Instan¬ 
taneous  field  amplitude  and  phase  front  angle  fluctuations.  The  valid¬ 
ity  of  these  assumptions  is  to  be  verified  a  posteriori. 


i.e.  one  equation  for  the  forward  coherent  amplitude  A,  one  for  the 
forward-scattered  irradiance  T*,  and  one  for  the  backscattered 
irradiance  I”.  The  coordinate  z  is  the  distance  along  the  beam  axis,  r 
is  the  radial  coordinate,  F  is  the  focal  distance  and  v|  is  the 
Laplacian  operator  in  the  plane  normal  to  the  z-axis.  The  quantities 
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a  ,  a  ,  a+,  a“,  D+  and  D~  are  respectively  the  molecular  and  aerosol 
mass 

absorption  coefficients,  the  aerosol  forward-  and  backscattering 
coefficients,  and  the  forward  and  backward  lateral  diffusion 
coefficients. 

A  simplifying  approximation  is  implicit  in  eqs.  4-6.  The  par¬ 
ticulate  medium  was  assumed  locally  homogeneous,  i.e.  homogeneous  over 
a  volume  of  the  size  of  the  beam  diameter.  This  implies  that  the  coef 
ficients  a's  and  D's  are  functions  of  the  coordinate  z  only. 

3.0  PROPAGATION  COEFFICIENTS 

The  parameters  a's  and  D's  of  eqs.  4-6  constitute  the  bulk  opti 
cal  scattering  properties  of  the  random  medium.  They  are  formally 
related  to  the  properties  of  the  particles  and  surrounding  medium 
through  the  closure  relations  derived  in  Ref.  15.  Those  closure  rela¬ 
tions  are  not  coupled  with  the  propagation  eqs.  4-6.  Hence,  the  prob¬ 
lem  of  determining  the  a's  and  D's  can  be  dealt  with  independently. 
However,  the  mathematics  is  difficult,  it  requires  solving  for  the 
covariance  function  of  the  geometrical-optics  phase  front  angle  fluctu 
atlons.  That  function  obeys  a  three-dimensional  nonlinear  partial 
differential  equation  with  coefficient  functions  of  disparate  scales. 
This  constitutes  a  formidable  problem  which  has  not  been  resolved  yet. 

Since  the  coefficients  a's  and  D's  are  independent  of  the  propa 
gation  equations,  it  appears  justified  to  seek  a  substitute  and  mathe¬ 
matically  more  convenient  method  to  determine  them.  The  most  suitable 
candidate  is  the  Mie  theory.  We  therefore  propose  in  this  report  to 
define  the  propagation  coefficients  heuristically  in  terms  of  conven¬ 
tional  statistical  summations  of  Mie  solutions  for  single-particle 
interactions.  The  method  is  based  on  the  physical  definitions  of  the 
coefficients  a’s  and  D's. 
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From  simple  examination  of  eq.  4,  we  find  that  the  coefficients 

a  ,  a  ,  a+  and  a~  model  the  fading  of  the  incident  coherent  beam 
m’  a*  s  s 

through  losses  caused  by  molecular  absorption,  aerosol  absorption, 
forward-  and  backscattering.  The  molecular  absorption  is  assumed  known 
or  given.  As  for  the  aerosol  terms,  the  descriptions  just  given  are 
consistent  with  the  standard  definitions  of  aerosol  absorption  and 
scattering.  Hence,  for  spherical  particles,  we  set 


x/2 


2n 


dQ_ 


i+  =  /  d9  sine  J  dip  J  dp  d^p)  xp2  ^-1  (9,H>,p), 


o  o 


x/2 


o  o 


3p" 

N(f 


[?] 


u  2x  ®  v  dQ 

a“  “  /  d0  sin0  f  dfr  /  dp  ■  up2  (e,4>,p),  [8] 


/  dP  "-app)  ’tP2  Qe<P> 


-  a+  -  a“ 


[9] 


where  0  is  the  scattering  angle,  4>  is  the  azimuthal  angle,  p  is  the 
particle  radius,  dN(p)/dp  is  the  particle  size  probability  density, 
dQ  (0,4>,p)/dQ  and  Q  (p)  are  the  differential  scattering  and  total 
extinction  efficiencies  of  a  particle  of  size  p,  and  dQ  «*  sin0dOd<t>  is 
the  differential  solid  angle.  The  size  probability  density  function  is 
a  specified  input.  The  scattering  and  extinction  efficiencies  Qg  and 
Qe  are  functions  of  the  radiation  wavelength,  particle  size,  and  com¬ 
plex  refractive  indices  of  the  particle  and  surrounding  medium.  They 
are  standard  expressions  of  the  Mie  solutions.  Hence,  the  coefficients 

a  ,  a+  and  a~  are  obtainable  through  a  series  of  operations  well  docu- 
a  s  s 

mented  in  the  literature. 


The  diffusion  coefficients  D+  and  D“  are  modeled  as  follows: 
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2n  n/2 

j  d4>  /  d9  sinQ  p(9,4>) 

D+(z-z’)  “  (z-z')  — - 2 -  ,  [lO] 

2n  n/2 

/  d<t>  /  de  p(e,$) 

o  o 

2n  n 

/  d$  /  d 0  sin9  p(9,4>) 

D“(zf-z)  -  (z'-z)  — - 2^1 -  ,  [11  ] 

2n  p(n) 

where  |z-z'|  is  the  distance  from  the  point  of  scattering  and  p(9,<J>)  is 
the  phase  function  defined  by 

p(9,<t>)  - - - -  /  dp  np2  (9,4>,p).  [12] 

(“s  +  ag)  o  dp 

The  linear  z-dependence  is  suggested  by  the  asymptotic  form  of  the 
closure  relations  of  Ref.  15  in  the  limit  of  small  jz-z'|.  The  propor¬ 
tionality  functions  are  averages  of  the  sine  of  the  scattering  angle  in 
the  forward  and  backward  hemispheres.  The  different  normalization 
constants  were  chosen  for  best  fit  with  data.  Although  the  formal 
closure  relations  indicate  that  D+  and  D~  should  saturate  with  increas¬ 
ing  | z-z’ | ,  calculations  have  shown  that  the  saturation  has  only  mar¬ 
ginal  effects  for  the  conditions  investigated  to  date. 

Equations  7-11  provide  a  model  for  calculating  the  coefficients 
of  the  propagation  eqs.  4-6.  These  formulas  were  obtained  by  applying 
the  results  of  the  Mie  theory  to  the  physical  processes  represented  by 
the  a's  and  D's.  The  connection  is  heuristic  but  analogous  to  what  is 
done  in  the  radiative  transfer  model.  Indeed,  the  radiative  transfer 
equation  is  a  mathematical  representation  of  the  scattering  phenomenon 
derived  from  the  conservation  of  flux.  It  involves  a  scattering  (or 
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phase)  function  which,  in  practice,  is  calculated  from  independent  Mie 
solutions  as  are  the  a's  and  D's  of  the  present  model. 

4.0  SOLUTIONS 

Equations  4-6  constitute  a  closed  system  for  the  solution  of  the 
forward  and  backscattered  irradiances  resulting  from  the  interaction  of 
a  coherent  beam  with  a  discrete  random  medium.  The  system  includes  an 
equation  for  the  amplitude  A  which  yields  the  reduced  irradiance  AA* 
and  two  coupled  equations  for  the  scattered  irradiance  functions  I"1"  and 
I-.  These  differential  equations  express  the  rates  of  change  of  A,  T* 
and  I-  due  to  scattering  losses  and  gains,  lateral  diffusion  and  dif¬ 
fraction.  The  scattering  losses  are  governed  by  the  extinction  coef¬ 
ficients.  For  1+  and  I“,  the  extinction  coefficient  (a  +  a  +  a-) 
lacks  the  term  a*  compared  with  that  for  the  reduced  Irradiance.  This 
means  that  the  forward-scattering  process  does  not  constitute  a  loss  to 
the  scattered  components  of  the  propagating  waves,  as  expected.  The 
scattering  gains  or  multiple-scattering  effects  are  given  by  the  right- 
hand  side  input  terms.  The  reduced  irradiance  feeds  Into  both  the 
forward-  and  backscattered  Irradiances  while  the  counterpropagating 
scattered  waves  feed  into  one  another.  Finally,  the  lateral  diffrac¬ 
tion  and  diffusion  broadenings  are  modeled  by  the  Laplacian  terms  of 
eqs.  4-6.  Diffraction  also  affects  the  broadening  of  I*  and  I",  but  it 
has  been  left  out  of  eqs.  5  and  6  as  beam  diffraction  is  much  smaller 
than  lateral  diffusion  by  scattering. 

The  only  difficulty  in  solving  the  system  of  eqs.  4-6  is  the 
coupling  between  the  forward-  and  backscattered  irradiances.  Here,  we 
proceed  Iteratively  by  first  setting  1“  ■  0  in  eq.  5  and  solving  for  I+ 
and  I-  in  sequence.  Higher  order  iterations  are  then  obtained  by 
substituting  the  I”  solution  back  into  eq.  5  and  repeating  the  process. 
As  it  turns  out,  the  first  Iteration  is  sufficient  in  many  practical 
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situations.  What  is  neglected  is  the  backscattering  of  the 
backseat tering. 

The  iterative  solution  can  easily  be  obtained  in  analytic  form 
if  the  original  beam  has  a  Gaussian  profile.  In  this  case,  the  first- 
iteration  solution  is 


AA* ( z , r ) 


P0  T(2) 


exp[- 


w2W(z) 

o 


]. 


[13] 


I+Cz.r) 


2 

/  dz' 
o 


F(z,2* ) 
H( z,z  ’ ) 


exp  [- 


w2H(z,z’) 


■]. 


P 

I-(z,r,Z)  «  — 

XW2 

o 


Z 


/  dz ' 


G(z, Z ' ) 
K( 2, 2 ' ) 


exp  [- 


W2K(z,2’) 


[14] 


P  Z  z’ 

+  °.  j  dz'  /  dz 

XWZ  z  o 

o 


L(z,z* ,z") 
M( Z , Z ' ,2*) 


exp 


[- 


w2M(z,z' ,z") 


•].[15] 


where  P  is  the  beam  power,  w  is  the  beam  radius  (1/e  in  irradiance) 
o  o 

at  the  boundary  of  the  cloud  (z  *0),  Z  is  the  depth  of  the  cloud,  and 
z  and  r  are  respectively  the  longitudinal  and  transversal  coordinates 
of  the  point  of  observation.  The  functions  T,  F,  G  and  L  are  "trans¬ 
mittance"  functions  and  W,  H,  K  and  M,  "beam-spreading"  functions. 

They  are  given  by  integrals  over  the  z-coordinate  of  the  propagation 
coefficients,  i.e. 


T(z)  -  exp{-/  [a  (z*)  +  a  (z')  +  a+(z')  +  a“(z,)]dz'},  [16] 

Ul  a  8  9 

o 
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F(z,z')  -  a+(z')  T(z’)  exp[-J  [a  (z")  +  a  (z”)  +  a“(z")]dz"},  [17 ] 

S  |  ul  B  S 

z 

z* 

G(z,z')  -  a“(z')  T(z')  exp{-J  [a  (z”)  +  a  (z")  +  a“(z") ]dz“ } ,  [18  ] 

5  mas 

z 

z' 

L(z,z’,z")  -  a-(z’)a+(z")T(z")  exp{-/  [b(z”  )  +  a  (z“*  )  +  a"(z~  ) ]dz”, 

S  S  ••  ul  a  S  i 

Z  * 

z' 

x  exp{-/  [am(z")  +  a^(z")  +  a~(z")]dz"},  [ 1 9 ] 

[20] 


w(z)  -  illill  +  _il_  , 

F2  k2w4 
o 


H(z,z’)  -  W(z’)  +—  J  <Z'~F>2  D+(z"-z»)dz", 

w2  z'  (z"-F)2 
o 

K(z.z')  -  W(z»)  +—  /  iilllli  D“(z’-z”)dz", 
w2  z  (z"-F)2 


M(z,z'  ,z")  -  W(z")  +  —  /Z  (Z"~F)2  D+(z-  -z”)dz~ 

v2  z"  (z"-F)2 

+  —j  (z'~F-)2  D-(zt-z")dz". 

w2  z  (z"-F) 2 
o 


[21] 


[22] 


[23] 


Equations  13-15  give  the  profiles  of  the  reduced  coherent  and 
forward-  and  backscattered  Incoherent  lrradiances  resulting  from  a 
Gaussian  beam  directed  into  an  inhomogeneous  particulate  medium.  These 
solutions  account  for  multiple  scatterings  of  all  orders. 
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5.0  FIELD-OF-VIEW  EFFECT 


The  solutions  given  by  eqs.  13-15  are  irradiance  solutions  and, 
therefore,  applicable  to  open  detectors  only.  For  detection  at  limited 
fields  of  view,  the  angular  distribution  of  the  radiation  flux  must  be 
taken  into  account.  Here,  we  model  this  effect  by  multiplying  the 
differential  contributions  to  the  scattered  radiation  in  eqs.  14  and  15 
by  a  correction  factor,  as  follows: 


F+Xz.r.Q^ 


P  z 
o 


w*  o 
o 


o  ’  w2H(z,z ' ) 


F*(z,r,Z,Q  )  - 

O  9 

nwz 
o 


P  Z 
o 


/  <!*•»<*•  ...r .0  >  «xp  [-  _£ - _] 

z  V  »  /  w2K(z,z  ) 


Z  z 

/  dz'v(z' tz,rfQ0)  /  dx-  .xp  [- 


w2M(z,z' 

o 


,z") 


-].  [25] 


where  F+(z,r,QQ)  and  F”(z,r,Z,Qo)  are  the  forward-  and  backscattered 
powers  detected  within  a  solid  angle  Qq  at  point  (z,  r),  and 
U(z',z,r,Q^  and  V(z',z,r,Qo)  are  f ield-of-view  reduction  functions  to 
be  specified  belcw.  Equation  13  does  not  need  modification  since  the 
coherent  component  conserves  its  small-divergence  characteristics. 


The  functions  U  and  V  are  reduction  factors  which  characterize  a 
receiver  of  field  of  view  Qq  relative  to  an  open  receiver.  In  other 
words,  U  and  V  in  their  respective  forward-  and  backscatter  configura¬ 
tion  measure  the  fraction  of  radiation  originating  from  a  scattering 
plane  z*  and  collected  by  a  receiver  of  field  of  view  compared  to 
that  collected  by  an  open  receiver.  U  and  V  are  therefore  smaller  than 
or  equal  to  unity.  In  the  limit  U«1  or  V*l,  eq.  24  or  25  becomes 
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identical,  by  definition,  to  the  original  irradiance  solution  14  or  15. 
With  reference  to  Fig.  1  for  a  schematic  definition  of  the  variables. 


U(z',z,r,Qo)  is  given  by 


U(z’,z,r,Qo)  -  e 


-  /  a+dzl"  /  d2r'  AA*(z’,r')  p(z’,0  )  R(0  ,Qq) 
z* 

/  d2r*  AA*(z',r*)  p(z',e  ) 


-  /  a+d2 


]  I  1- 

m*l 


fdz  P(m  z  ) 
m  1 


JdV  /d2r1.../d2rraR(er,Qo)AA*(z’,r,)p(z,,eo)^niP(z1,ei) 


Jd2r,fd2Tl...Jd2rm  AA*(z' ,r,)p(z,,0o)  I^pCz^Gj) 


.  [26] 


where  p  is  the  phase  function,  R  is  the  receiver  function,  the  0's  are 


scattering  angles  which  can  be  expressed  in  terms  of  z,  z' ,  z^ ,  r'  and 


r^  by  simple  trigonometry,  PCmJz^)  is  the  joint  probability  of  a 


particular  configuration  z^  of  m  scattering  planes,  and  M  is  the  maxi¬ 


mum  number  of  scattering  planes.  A  similar  expression  for  V(z',  z,  r. 


Qq)  is  obtained  by  replacing  AA*  by  I*  in  eq.  26. 


Equation  26  for  U  has  two  main  contributions.  The  first  term  is 
the  contribution  from  the  radiation  that  is  left  unscattered  following 
the  primary  scattering  at  z'f  and  the  second  term  is  the  contribution 
from  the  radiation  which  has  experienced  at  least  one  scattering 
between  z'  and  z.  For  both  terms  the  fraction  of  detected  power  is 
calculated  by  integrating  over  all  scattering  possibilities  and 


ratioing  the  result  for  a  receiver  function  R(0r,QQ)  to  that  for  an 


open  receiver,  i.e.  R«l. 


FIGURE  1  -  Schematic  representation  of  scattering  events  for 

calculating  the  f ield-of-vlew  reduction  factors  U  and 
a)  single  scattering;  b)  multiple  scatterings 
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The  first  term  of  eq.  26  is  straightforward.  The  numerator 
integral  is  proportional  to  the  quantity  of  radiation  scattered  from 
plane  z'  which  falls  within  a  field  of  view  of  Q  at  (z,r)  after  propa- 
gating  the  distance  (z-z*)  without  scatterings.  The  denominator  is  the 
corresponding  result  for  an  open  detector.  The  second  term  is  more 
complex  but  the  principle  remains  the  same.  The  surface  integrals  give 
a  quantity  proportional  to  the  amount  of  radiation  detected  at  (z,r) 
after  being  scattered  m  times  between  z'  and  z,  for  a  receiver  function 
R(9r,Qo)  in  the  numerator  and  for  an  open  receiver  in  the  denominator. 
The  contributions  are  summed  over  all  orders  of  scattering.  As  indi¬ 
cated  by  the  expressions  in  the  numerator  and  denominator,  the  differ¬ 
ential  scattered  radiation  reaching  the  detector  is  proportional  to  the 
m 

product  p(z',e  )  n  p(z.,0.).  At  optical  and  infrared  wavelengths  and 
°  i-1  1  1 

for  atmospheric  aerosols,  the  phase  functions  pCz^,©^)  are  peaked  in 
the  forward  direction.  Hence,  the  resulting  integrals  decrease  rapidly 
with  i  for  i  >  1  unless  the  0^’s  are  near  zero.  It  follcws  that  the 
ratio  of  the  surface  integrals  in  the  last  term  of  eq.  26  is  about  the 
same  for  all  values  of  m.  In  other  words,  the  m  scatterings  can  be 
approximated  for  the  purpose  of  calculating  this  ratio  by  a  single 
scattering  with  non-negligible  angle.  The  ratio  can  thus  be  factored 
out  of  the  summation  and  since  that  one  scattering  can  occur  with  equal 


probability  anywhere  between  z*  and  z,  we  approximate 


I'l  4L»  Vi  | 

I  J  dz|  . . .  / dz  P(m|z  )  R(z.)  *  /  - RCzj), 

»-l  “  1  1  z’  z-z’ 


where  R  denotes  the  ratio  of  the  surface  integrals.  We  thus  obtain 
from  eq.  26  the  following  practical  expression  to  calculate  the  field- 


of-view  reduction  factor  U: 


U(z',z,r,Qo)  *  e 


-  J  a+dz  /dVAA*(z',r')p(z\eg)R(er,Qo) 


Jd2r,AA*(z’ ,r')p(z'  ,0  ) 


y  / 
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/ d2r'  Jd2£]R(0r  ,Qo)AA*(z*  ,r' )p(z'  ,eQ)p(z1 ,0!  ) 


(z-z')  z*  |^/d2r' Jd2r1AA*(z' ,r')p(z' ,0o)p(z1,e1) 

The  corresponding  factor  V(z',z,r,2o)  is  obtained  by  substituting  I+ 
for  AA*. 


Equation  28  completes  the  method  of  solution  for  the  forward  and 
backward  flux  density  profiles  (AA*  +  F*)  and  F_  for  a  given  receiver 
geometry.  The  required  inputs  are  the  receiver  function  R(0r,QQ)  and 
the  Mie-derived  phase  function  p  and  absorption,  scattering  and  diffu¬ 
sion  coefficients  a's  and  D's.  The  reduction  factors  U  and  V  are 
calculated  by  substitution  in  eq.  28  of  the  irradiance  solutions  13  and 
14  for  AA*  and  I*.  Finally,  the  required  AA*,  F*  and  F“  are  obtained 
from  eqs.  13,  24  and  25. 


For  the  general  case,  the  determination  of  the  reduction  factors 
U  and  V  is  computation-intensive.  However,  the  surface  integrals  of 
eq.  28  can  be  carried  out  analytically  at  considerable  savings  of  com¬ 
putation  effort  if  AA*,  I*,  p  and  R  are  approximated  by  Gaussian  func¬ 
tions.  Equation  13  shows  that  this  is  exactly  the  case  for  AA*  while 
eq.  14  gives  I*  as  the  sum  of  Gaussian  profiles,  which  satisfies  the 
condition  for  analytic  integration.  The  phase  function  is  not  Gaussian 
but  the  forward  peak,  which  is  the  dominant  contribution  here,  can  be 
modeled  by  a  Gaussian  function  with  reasonable  accuracy.  In  any  case, 
it  is  the  width  of  the  peak  that  really  matters  in  the  evaluation  of  U 
and  V  and  not  the  exact  functional  form.  Finally,  the  function 
R(0r,Qo)  °f  oany  practical  receivers  is  well  represented  by  a  Gaussian. 
We  therefore  use  the  Gaussian  approximation  to  perform  the  surface 
Integrals  of  eq.  28  and  it  is  hoped  that  the  ratloing  operation  will 
smooth  out  the  differences. 


mm 


UNCLASSIFIED 

16 


If  we  let 


AA*(z’,r')  or  I+(z',r')  ~  e 


-  r*  /w2(z') 


-  e2/o2Cz') 


p(z’,9g)  ~  e 


p(zi,e-)  ~  e 


-  6$/cj2(zi) 


R(e  ,Q  )  ~  e 

r*  o 


-e2/w2 

r  o 


and  if  we  make  the  paraxial  approximation,  i.e.  the  0's  <  0.2,  we 
obtain 

/d2r'AA*(z' ,r,)p(z’,es)R(er,Q0)  ^  +W2(z»)b-J 

fd2r,AA*(z,,r’)p(z,,0  )  Ll  +  w2(z ' )a J 

X  exp  j (" - - - -  - - - lr2j, 

||_1  +  w2(z')b  1  +  w2(z')aj  i 

/d2r'/d2r1AA*(z' ,r')p(z' ,98)p(zi,©i)R(9r,C0)  rd 


AA*(z' 
c2f  2 


(dg  -  f 2 )(de 


r*  i 


iCW 


Equations  33  and  34  are  substituted  in  eq.  28  to  calculate  the  reduc¬ 
tion  factor  U  or  V.  The  function  w2  is  determined  from  the  irradiance 
solution  13  or  14,  o2  from  the  phase  function  and  a>2  is  assumed  to  be 
given.  What  remains  to  be  done  numerically  is  the  single  integration 
over  It  is  straightforward  to  verify  that  in  the  limit  u>2  >>  a2, 

i.e.  for  an  open  receiver,  the  factors  U  and  V  tend  to  unity  as 
required. 

The  model  given  by  eqs.  13,  24  and  25  for  the  flux  density 
profiles,  eqs.  7-11  for  the  absorption,  scattering  and  diffusion 
coefficients  a’s  and  D*s,  and  eqs.  28,  33  and  34  for  the  f ield-of-view 
reduction  factors  U  and  V  has  been  coded  for  numerical  applications. 
The  program  calculates  the  transmitted  and  backscattered  profiles  for 
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specified  detectors  and  propagation  media  which  can  be  inhomogeneous. 
The  code  requires  only  limited  memory  space  and  computation  time. 

6.0  COMPARISON  WITH  MEASUREMENTS 


6.1  Experiment 

For  comparison  with  data,  an  experiment  was  designed  to  measure 
transmission  through  and  backscatter  from  well-defined  clouds.  The 
most  convenient  way  of  achieving  reliable  control  on  cloud  characteris¬ 
tics  is  to  generate  them  in  an  enclosed  laboratory  environment  where 
the  parameters  can  be  varied,  monitored  and  maintained  constant.  The 
experiment  was  carried  out  by  Optech  Inc.  under  contract  for  DREV 
(Ref.  16).  Water  droplet  clouds  were  generated  by  ultrasonic  nebu¬ 
lizers  in  a  chamber  of  adjustable  length.  The  size  distribution  of  the 
drops  produced  by  those  nebulizers  was  reasonably  invariant.  The  cloud 
density  was  varied  by  regulating  the  nebulizers  output  flow  into  and 
out  of  the  chamber.  After  sufficient  settling  time  and  through  gentle 
stirring,  a  homogeneous  concentration  was  maintained  for  as  long  as 
needed. 


The  size  distribution  of  the  water  drops  was  measured  with  a 
Classical  Aerosol  Scattering  Probe  (CASP)  manufactured  by  Particle 
Measuring  Systems  Inc.  and  modified  for  high  concentration  measure¬ 
ments.  The  drop  diameters  were  distributed  between  approximately  0.1 
and  15  pm  with  a  peak  near  1.0  pm.  There  are  some  uncertainties 
regarding  the  accuracy  of  the  CASP  probe.  In  particular,  it  measured 
more  small  particles  than  expected  from  previous  experience.  However, 
as  no  alternate  method  of  measurement  could  be  explored  because  of 
resource  constraints,  and  since  the  results  were  reproducible  for  most 
cloud  densities,  the  measured  distributions  were  considered  representa¬ 
tive  of  the  true  size  distribution  of  the  nebulizer-generated  clouds. 

A  typical  distribution  is  illustrated  in  Fig.  2. 
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Measurements  were  performed  at  three  wavelengths:  0.63,  1.06 
and  10.6  pm.  The  extinction  coefficient  was  continuously  monitored  by 
recording  the  cloud  transmittance  at  0.63  or  1.06  pm.  Empirical 
relations  had  previously  been  established  between  extinction  at  the 
three  wavelengths  of  the  experiment.  Hence,  the  relevant  extinction 
coefficient  could  easily  be  determined  from  the  measured  0.63  or 
1.06-pm  value.  A  range  of  extinction  values  (at  0.63  pm)  between  0  and 
3m-1  was  covered  for  a  cloud  depth  of  1.5  or  2.6  m.  The  remaining 
parameters  necessary  for  model  predictions  were  determined  by  Mie 
calculations  based  on  the  size  distribution  shown  in  Fig.  2  and  the 
complex  refractive  index  of  water  at  0.63,  1.06  and  10.6  pm.  Table  I 
summarizes  the  relative  parameter  values  for  the  conditions  of  the 
experiment. 

6.2  Transmitted  Beam  Profiles 


The  transmitted  beam  profiles  were  obtained  at  the  three 
wavelengths  of  the  experiment.  They  were  measured  by  scanning  a 
receiver  of  variable  field  of  view  a  short  distance  from  the  exit 
window  of  the  cloud  chamber.  Optical  attenuation  was  used  to  increase 
the  dynamic  range  of  the  detectors,  and  spurious  side-lobe  variations 
were  removed  by  spatially  filtering  the  beam.  The  unscattered  profiles 
at  0.63,  1.06  and  10.6  pm  are  shown  in  Figs.  3-5.  As  illustrated,  the 
profiles  are  well  approximated  by  a  Gaussian  function  over  about  4 
orders  of  magnitude.  Beam  diameters  at  1/e2  in  irradiance  are 
respectively  1.15,  0.90  and  1.60  cm  for  the  0.63,  1.06  and  10.6-pm 
beams . 


The  profiles  were  measured  for  various  extinction  coefficients 
at  two  fields  of  view,  20  and  350  mrad,  and  for  two  cloud  depths,  1.5 
and  2.6  m.  Figure  6  compares  data  with  model  predictions  for  the 
0.63-pm  beam,  20  and  350-mrad  fields  of  view,  1.5-m  cloud  depth,  and 
extinction  coefficient  of  about  3.  The  agreement  is  excellent.  The 


W 
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central  core  of  the  beam  keeps  its  Gaussian  shape  and  appears  unaffect¬ 
ed  by  multiple  scatterings.  This  is  the  case  for  both  fields  of  view 
and  it  is  well  confirmed  by  the  model  predictions  which  reveal  that  the 
central  core  consists  almost  exclusively  of  the  reduced  irradiance. 

Away  from  the  center,  the  profiles  suddenly  level  off  to  show  the  beam¬ 
broadening  effect  caused  by  multiple  scatterings.  As  illustrated  in 
Fig.  6,  the  broadening  is  a  function  of  the  receiver  field  of  view. 
There  is  a  difference  of  more  than  2  orders  of  magnitude  between  the 
off-axis  levels  recorded  at  20  and  350  mrad.  The  broadened  profiles  at 
both  20  and  350  mrad  agree  very  well  with  the  calculated  curves.  The 
corresponding  results  at  1.06  pm  are  plotted  in  Fig.  7.  The  agreement 
for  the  350-mrad  data  is  not  as  good  as  in  Fig.  6  but  is  still  within 
expected  experimental  errors. 

The  1.06-pm  profile  measurements  were  repeated  for  a  cloud  depth 

of  2.6  m.  The  results  are  shown  in  Fig.  8.  The  agreement  at  350  mrad 

Is  good  despite  some  asymmetry  in  the  measured  profile,  but  the  calcu¬ 
lated  off-axis  level  at  20  mrad  is  greater  than  the  experimental  values 
by  a  factor  of  about  5.  It  is  difficult  to  explain  the  origin  of  this 
discrepancy,  but  a  slight  detector-axis  misalignment  of  the  order  of  10 
mrad  with  respect  to  the  beam  axis  could  cause  a  signal  drop  of  this 
magnitude. 

At  10.6  pm.  It  was  very  difficult  to  observe  any  beam  broadening 

as  the  forward  scattering  peak  is  much  wider  because  of  the  long  wave¬ 

length  compared  with  the  particle  radii.  The  only  conditions  that 
revealed  measurable  off-axis  signal  were  for  a  field  of  view  of  350 
mrad  and  a  cloud  depth  of  2.6  m.  These  data  are  plotted  in  Fig.  9 
where  it  is  shown  that  the  agreement  with  the  calculated  solution  is 
nearly  exact. 
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Figures  10-12  give  the  transmitted  profiles  measured  with  a 
350-rarad  field  of  view  for  various  extinction  coefficients  between  0 
and  3.  The  model  solutions  are  well  confirmed  in  the  central  core  for 
all  conditions  illustrated.  The  solutions  of  Figs.  10  and  11  show  that 
the  off-axis  multiscattering  level  and  drop  rate  are  weak  functions  of 
the  extinction  coefficient  for  the  1.5-ra  cloud  in  the  range  investigat¬ 
ed.  This  is  well  verified  in  Fig.  11  for  the  1.06-pm  beam  although  the 
measured  values  are  higher  than  the  model  predictions  by  a  factor  of 
about  3.  In  Fig.  10  for  the  0.63-pm  beam,  there  is  a  non-negligible 
offset  between  the  off-axis  measurements  at  0.8  and  1.8  m-1  and  those 
at  3.2  m-1.  This  offset  is  not  borne  out  by  the  model  but  may  be 
explained  by  experimental  errors  as  the  data  points  show  a  rather  large 
asymmetry  in  that  case.  Finally,  the  solutions  of  Fig.  12  for  a  cloud 
depth  of  2.6  m  indicate  a  greater  variation  of  the  off-axis  level  with 
extinction  coefficient  in  agreement  with  the  data. 

In  summary,  the  transmitted  beam  profile  solutions  of  the  propo¬ 
sed  model  are  well  corroborated  by  the  laboratory  measurements  of  Ref. 
16.  The  parameters  were  varied  over  a  wide  range:  three  wavelengths 
from  the  visible  at  0.63  pm  to  the  infrared  at  10.6  pm;  extinction 
coefficients  between  0  and  3  m-i  for  a  cloud  depth  of  1.5  or  2.6  m 
which  corresponds  to  optical  depths  between  0  and  7;  and  two  widely 
separated  receiver  fields  of  view  of  20  and  350  mrad.  In  all  cases, 
the  predicted  profiles  reveal  a  central  core,  which  retains  its 
Gaussian  shape  and  is  essentially  the  reduced  coherent  irradiance, 
flanked  by  a  slowly  decreasing  flux  density  level  due  to  forward  multi¬ 
ple  scatterings.  The  measured  magnitude  of  this  broadened  flux  density 
depends  strongly  on  the  receiver  field  of  view.  Those  general  features 
of  the  solutions  are  in  excellent  agreement  with  the  data.  There  are 
some  discrepancies  on  the  off-axis  levels  but  it  cannot  be  ascertained 
if  they  result  from  experimental  or  model  errors.  In  any  case,  there 
is  no  trend  nor  pattern  that  could  point  to  a  systematic  flaw  in  the 
proposed  model. 


NORMALIZED  POWER  NORMALIZED  POWER 
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6.3  Backscatter 


The  lidar  technique  has  great  potential  as  a  practical  remote 
sensing  device  for  determining  the  extinction  coefficient  of  suspended 
natural  and  artificial  aerosols.  However,  most  existing  inversion 
methods  neglect  the  multiple-scattering  contributions,  which  may  have 
drastic  adverse  effects  on  the  outcome.  The  present  backscatter  solu¬ 
tion  (eq.  25)  could  help  resolve  that  problem.  Indeed,  since  eq.  25  is 
written  in  analytic  form,  it  can  in  principle  be  inverted  to  retrieve 
the  extinction  coefficient  from  lidar  returns  affected  by  multiple 
scatterings.  It  is  therefore  important  to  validate  eq.  25.  For  that 
reason,  the  Optech  experiment  (Ref.  16)  was  also  designed  to  generate 
relevant  backscatter  data. 

The  backscatter  receiver  consisted  of  a  spherical  mirror  aper- 
tured  to  a  5-cm  diameter.  The  field  of  view  was  controlled  by  a  field 
stop  placed  in  the  image  plane  corresponding  to  the  cloud  chamber. 
Fields  of  view  of  10  and  15  mrad  were  chosen;  15  mrad  is  the  maximum 
safe  value,  given  the  size  of  the  chamber  windows  and  cloud-to-receiver 
separation  distance.  A  dichroic  beam  splitter  mounted  behind  the  field 
stop  allowed  simultaneous  detection  at  10.6  pm  and  either  0.63  or  1.06 
pm.  The  receiver  was  positioned  as  close  as  possible  to  the  laser  axis 
and  aligned  for  total  beam  overlap  over  the  length  of  the  cloud.  The 
main  experiment  was  conducted  for  a  cloud-to-receiver  distance  of  10  m 
but  additional  tests  were  made  at  8.8  and  5.9  m. 

The  backscatter  returns  from  the  integrated  cloud  length  were 
measured  as  functions  of  wavelength,  extinction  coefficient,  receiver 
field  of  view,  cloud  depth  and  cloud-to-receiver  distance.  The  cloud 
return  signals  were  normalized  to  the  signal  reflected  from  a  reference 
Lambertian  target  placed  in  front  of  the  chamber  window  at  the  begin¬ 
ning  of  each  measurement  period.  A  barium  sulfate  (BaSOt,)  target 
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was  used  at  0.63  and  1.06  pm  and  a  sanded  styrofoam  panel  at  10.6  pm. 

A  reflectivity  of  100%  (Ref.  17)  is  assumed  for  the  barium  sulfate,  and 
35%  (Ref.  18)  for  the  styrofoam.  The  dimensional  backscatter  can  then 
be  calculated  by  multiplying  the  normalized  data  by  PQr2v/L2  where  Pq 
is  the  laser  power,  r^  the  receiver  radius,  v  the  target  reflectivity 
and  L  the  target-to-receiver  distance.  However,  the  reported  normal¬ 
ized  backscatter  results  at  0.63  and  1.06  pm  are  smaller  than  1.06-pm 
data  obtained  earlier  under  similar  conditions  (Ref.  19).  A  systematic 
error  in  the  measurement  of  the  reference  target  return  appears  to  be 
the  cause.  Subsequent  tests  (Ref.  20)  performed  with  two  different 
receivers  confirmed  that  the  measured  values  of  Ref.  16  are  indeed  too 
small  by  a  multiplicative  factor  between  2  and  3.  To  settle  the  prob¬ 
lem,  we  decided  to  determine  the  calibration  constant  by  fitting  the 
measurements  in  the  small-extinction  limit  to  calculations  based  on  the 
integration  of  the  single-scattering  lidar  equation,  i.e. 


P(Z)  *  Pq - P(n,Z)  exp  {-2  /  ag(z)dz}, 

Z2  o 


[42] 


where  P(Z)  is  the  range-resolved  lidar  return  per  unit  length,  Z  is  the 
range,  Pq  is  the  transmitter  power,  S  is  the  area  of  the  receiver 
aperture,  P(x,Z)  is  the  backscatter  coefficient,  and  ae(z)  is  the 
extinction  coefficient.  This  method  is  acceptable  since  the  single¬ 
scattering  lidar  equation  is  independent  of  our  propagation  model.  It 
has  yielded  a  correction  factor  of  3.15  for  both  the  0.63  and  1.06-pm 
data,  which  is  in  reasonable  agreement  with  the  correction  estimated 
from  the  post-experiment  tests. 


The  backscatter  measurements  at  0.63  pm  for  the  1.5-m  cloud, 

10-m  cloud-to-receiver  distance  and  10  and  15-mrad  fields  of  view  are 
compared  with  the  model  predictions  in  Fig.  13.  The  agreement  is 
almost  exact  except  for  the  15-mrad  case  where  the  theoretical  solution 
is  about  5%  higher  than  the  data  at  high  extinction  values.  These 
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small  differences  are  well  within  the  experimental  uncertainties  espe¬ 
cially  since  the  field-of-vlew  functions  were  not  directly  measured  but 
inferred  from  the  field-stop  sizes.  The  solution  obtained  by  integrat¬ 
ing  the  single-scattering  lidar  equation  is  also  drawn  for  comparison. 
It  is  in  perfect  agreement  with  the  present  solutions  for  extinction 
values  smaller  than  about  0.3  m_1  at  which  point  it  begins  to  saturate 
while  the  measured  returns  and  multiple-scattering  solutions  continue 
to  grow,  although  at  a  gradually  slower  rate.  These  results  show  that 
the  multiple-scattering  effects  indicated  by  the  f ield-of-view  depend¬ 
ence  and  larger  backscatter  signal  are  very  important  for  the  condi¬ 
tions  of  Fig.  13  and  well  accounted  for  by  the  present  model. 

The  corresponding  results  obtained  at  1.06  pm  are  plotted  in 
Fig.  14.  The  fit  is  not  as  close  as  in  Fig.  13  but  still  within 
acceptable  limits.  The  multiple-scattering  effects  are  less  than  those 
at  0.63  pm,  which  is  in  agreement  with  the  predictions.  In  particular, 
the  measured  f ield-of-view  differences  are  not  as  explicit  as  in  Fig. 
13,  but  it  must  be  noted  that  there  is  only  one  set  of  data  at  15  mrad 
for  three  independent  sets  at  10  mrad.  Thus,  the  10-mrad  experimental 
scatter  resulting  from  measurements  on  different  occasions  may  be 
responsible  for  confusing  the  f ield-of-view  discrimination.  There  is  a 
slight  change  in  the  data  slope  at  an  extinction  value  of  about  1.5 
m”1.  It  could  not  be  determined  if  this  change  is  real  or  is  an  arti¬ 
fact  of  the  experiment. 

The  results  obtained  at  10.6-pm  are  plotted  in  Fig.  15.  As  it 
turns  out,  the  water  clouds  are  much  less  efficient  scatterers  at  10.6 
pm  than  at  0.63  or  1.06  pm;  the  returns  are  smaller  by  about  2  orders 
of  magnitude.  This  explains  the  larger  experimental  scatter  of  the 
experimental  values  of  Fig.  15  compared  with  those  of  Fig.  13  or  14. 
Within  that  scatter,  the  low-extinction  half  of  the  data  agrees  very 
well  with  the  model  solutions.  However,  there  is  a  sudden  increase 
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near  the  extinction  value  of  0.4  m~ 1  which  is  not  borne  out  by  the 
solutions.  This  effect  is  probably  artificial  and  due  to  reflections 
of  forward-scattered  radiation  off  chamber  walls  (Ref.  20).  On  the 
other  hand,  there  is  very  good  agreement  on  the  weak  f ield-of-view 
dependence  between  the  10  and  15-mrad  receivers.  This  weak  dependence 
results  from  the  wide  angular  forward  scattering  demonstrated  by  the 
theoretical  and  experimental  results  given  in  Fig.  9.  In  other  words, 
there  are  very  few  multiple-scattering  contributions  within  fields  of 
view  of  the  order  of  10-15  mrad.  This  is  further  corroborated  by  the 
fact  (illustrated  in  Fig.  15)  that  the  single-scattering  solution  is 
only  slightly  less  than  the  present  solutions  for  10  and  15  mrad. 

The  cloud  depth  effect  was  also  investigated.  The  results  for  a 
cloud  depth  of  2.6  m  are  reproduced  in  Figs.  16  for  1.06  pm  and  17  for 
10.6  pm.  The  cloud-to-receiver  distance  was  reduced  to  8.8  m.  The 
measurements  at  1.06  pm  fit  the  solution  very  well  in  the  limit  of 
small  extinction  coefficients,  but  they  suddenly  change  slope  at  about 
0.4  m-1  to  fall  below  the  theoretical  curve  although  they  remain 
approximately  parallel  to  it  thereafter.  The  differences  are  of  the 
order  of  20-30%.  No  exact  explanation  could  be  found  but  it  is  most 
likely  an  experimental  error  since  the  high-extinction  returns  for  this 
case  are  nearly  equal  to  (or  even  smaller  than)  those  of  Fig.  14  for  a 
smaller  cloud  depth  and  greater  cloud-to-receiver  distance.  The  corre¬ 
sponding  results  at  10.6  pm  shown  in  Fig.  17  follow  approximately  the 
same  pattern  as  for  the  1.5-m  cloud  (Fig.  15).  There  is  a  fair  agree¬ 
ment  at  the  small-extinction  end,  but  the  calculations  are  smaller  than 
the  data  by  about  30%  following  the  sudden  increase  of  the  measurements 
near  0.4  m-1.  As  in  Figure  15,  the  same  explanation  of  contamination 
by  wall  reflections  is  probably  applicable. 

The  influence  of  the  cloud-to-receiver  distance  was  studied  by 
moving  the  2.6-m  chamber  closer  to  the  receiver,  from  8.8  to  5.9  m. 
Measurements  in  that  case  were  complicated  by  the  shorter  separation 


JU  JIUIUI'JIUI 


UNCLASSIFIED 

32 


and  longer  cloud  length,  which  made  it  more  difficult  to  adjust  for 
proper  overlap  and  caused  uncertainties  in  the  f ield-of-view  definition 
since  the  depth  of  focus  could  not  accommodate  the  whole  length  of  the 
cloud.  The  results  are  plotted  in  Figs.  18  for  1.06  pm  and  19  for 
10.6  pm.  In  both  cases,  the  agreement  is  good  at  small  extinction 
coefficients  but  the  subsequent  rise  of  backscatter  with  extinction 
coefficient  is  faster  than  predicted.  At  10.6  pm,  this  appears  to  be 
the  same  behavior  as  in  Fig.  17  for  8.8-m  distance  except  that  it  is 
amplified.  At  1.06  pm,  the  experimental  slope  is  now  greater  than  that 
of  the  theoretical  curve.  There  are  no  clear  explanations  for  these 
effects  but  a  probable  cause  may  be  related  to  the  overlap  and  depth- 
of-focus  problems  discussed  above.  If  indeed  present,  these  translate 
into  a  range-dependent  field  of  view  which  is  not  accounted  for  by  the 
fixed  function  used  for  the  calculations. 


In  summary,  the  backscatter  solutions  of  the  proposed  model  are 
well  validated  by  the  laboratory  data  of  Ref.  16.  The  basic  results 
given  in  Figs.  13-15  for  the  1.5-m  cloud  demonstrate  an  agreement  of 
better  than  10%  for  the  three  wavelengths  and  the  two  fields  of  view 
except  for  the  high-extinction  measurements  at  10.6  pm  which  are  prob¬ 
ably  contaminated  by  wall  reflections.  The  data  of  Figs.  13-15  are  the 
more  reliable  of  the  series  since  the  1.5-m  cloud/10-m  cloud-to- 
receiver  separation  constitutes  the  primary  configuration  of  the  Optech 
experiment  which  has  been  extensively  tested.  The  additional  data  with 
the  2.6-m  cloud  and  separation  distance  of  8.8  and  5.9  m  are  probably 
affected  by  depth-of-f ocus  and  overlap  problems.  There  are  obvious 
errors  such  as  the  smaller  returns  at  8.8  m  than  at  10.0  m  or  the 
larger-than-  expected  returns  for  high  extinction  coefficients  at  10.6 
pm.  Considering  these  effects,  the  theoretical  fits  of  Figs.  16-19 
appear  reasonable.  Moreover,  the  solutions  of  Figs.  16-19  exhibit 
features  easily  explained  by  simple  physical  arguments.  First,  there 
is  the  overall  increase  of  the  signal  level  with  a  decrease  of  the 
cloud-to-receiver  separation  which  is  simply  the  geometric-range 
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FIGURE  13  -  Measured  and  calculated  backscattered  power  from  a  water 
droplet  cloud  as  a  function  of  the  cloud  extinction 
coefficient  and  receiver  field  of  view.  The  curves  are 
calculations  and  the  symbols,  measurements.  SSCAT  is 
single-scattering  solution.  The  wavelength  is  0.63  pm;  the 
cloud  depth,  1.5  m;  and  the  cloud-to-recelver  separation. 
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factor.  Second,  there  Is  the  steeper  slope  at  small  extinction  coef¬ 
ficients  for  larger  cloud  depths.  This  occurs  because  at  small  extinc¬ 
tion  the  backscatter  from  the  entire  length  of  the  cloud,  even  for  an 
extended  depth,  contributes  to  the  detection;  whereas  at  high  extinc¬ 
tion  the  backscatter  from  the  far  end  is  more  attenuated  if  the  cloud 
is  longer.  Thus,  the  backscatter  is  proportional  to  cloud  length  in 
the  low-extinction  limit  but  gradually  becomes  independent  of  it  as  the 
extinction  grows.  This  explains,  in  agreement  with  the  present  model, 
why  the  integrated  backscatter  initially  increases  faster  for  greater 
cloud  depths. 
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FIGURE  19  -  Same  as  Fig.  18  except  that  the  wavelength  is  10.6  jim. 
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7.0  CONCLUSION 


The  proposed  multiple-scattering  propagation  model  is  well 
validated  by  the  laboratory  simulation  data  of  Ref.  16.  The  predicted 
flux  density  profiles  of  the  transmitted  radiation  show  a  central  core 
surrounded  by  a  low-level  and  much  more  uniform  flux  distribution. 

These  features  were  well  verified  over  the  complete  parameter  range 
investigated  which  included  three  wavelengths,  two  receiver  fields  of 
view,  and  optical  depths  between  0  and  approximately  7.  The  discrepan¬ 
cies  between  measured  and  predicted  off-axis  levels  show  no  particular 
trend  and  are  within  expected  measurement  errors. 

The  backscatter  solutions  are  in  excellent  agreement  with  the 
basic  1.5-m  cloud  data  for  three  wavelengths,  two  fields  of  view,  and 
optical  depths  between  0  and  A. 5.  For  the  larger  2.6-m  cloud  and 
shorter  cloud-to-receiver  separations  which  stretch  to  the  limit  the 
designed  receiver  specifications,  there  are  differences  of  20  to  100%. 
However,  considering  all  possible  causes  of  experimental  errors,  the 
fit  is  acceptable.  The  solutions  are  found  to  adjust  to  the  greater 
cloud  depth  and  shorter  separation  in  ways  consistent  with  simple  phys¬ 
ical  explanations,  in  particular,  the  more  rapid  rise  of  backscatter  at 
small  extinction  coefficients  and  the  overall  increase  in  the  signal 
level. 


A  complementary  test  was  performed  by  comparing  the  backscatter 
solutions  of  the  present  model  with  Monte  Carlo  calculations  (Ref.  21). 
The  computations  were  carried  out  for  a  1-km-deep  cloud  and  a  cloud-to- 
receiver  separation  of  200  m.  Three  wavelengths  (1.06,  3.5  and  10.5 
pm),  six  cloud  models  (maritime  aerosols  70  and  99%  RH,  rural  aerosols 
70  and  99%  RH,  and  urban  aerosols  70  and  99%  RH;  Ref.  22),  three 
extinction  coefficients  (0.5,  2.0  and  5.0  km-1),  and  two  fields  of  view 
(1  and  10°)  were  investigated.  For  each  case,  the  ratio  of  the  Monte 
Carlo  to  the  present  model  backscatter  returns  was  calculated.  The 
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resulting  average  ratio  is  1.01  with  a  standard  deviation  of  0.05,  and 
the  maximum  and  minimum  values  are  1.10  and  0.75  respectively.  There¬ 
fore,  the  agreement  between  the  two  different  solution  methods  is 
excellent  and  provides  further  corroboration  for  the  proposed  model  on 
a  scale  representative  of  atmospheric  applications. 

The  results  of  this  report  show  that  the  model  defined  by  the 
propagation  eqs.  4-6,  the  parameter  definitions  7-11,  and  the  field-of- 
view  functions  U  and  V  of  eq.  28  is  a  valid  representation  of  the 
multiple-scattering  effects  on  laser  propagation  and  scattering  in 
aerosol  clouds.  Forward-  and  backseat tering  phenomena  at  optical  and 
infrared  wavelengths  in  aerosols  of  extinction  coefficients  as  high  as 
3  m-1  were  properly  modeled.  This  constitutes  a  posteriori  confirma¬ 
tion  of  the  hypotheses  leading  to  the  propagation  eqs.  4-6  and  of  the 
heuristic  determination  of  the  scattering  and  diffusion  coefficients 
ct's  and  D’s.  The  solutions  are  given  by  analytic  expressions,  eqs. 
13-23,  which  handle  Inhomogeneous  clouds,  are  easily  coded,  and  require 
reasonable  computation  efforts.  Backscatter  calculations  are  faster 
than  with  the  biased  Monte  Carlo  code  of  Ref.  21  by  a  factor  of  at 
least  20.  Since  the  solutions  are  analytic,  they  are.  In  principle, 
applicable  to  the  inverse  lidar  problem  In  the  presence  of  multiple¬ 
scattering  contributions.  This  is  an  interesting  property  that  will  be 
investigated  In  future  work. 
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